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I discuss a method to model the instrumental response of the CMS and ATLAS detectors at high 
missing transverse energy to dominant standard model V+jets backgrounds, where 1/ is a 7 or 
W , using multi-jet QCD events. The method is developed for new physics searches in early data at 
the Large Hadron Collider (LHC) with minimal recourse to simulation. 
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I. INTRODUCTION 

The LHC enters a new energy regime to explore the 
origin of the electroweak symmetry breaking, search for 
and study physics beyond the standard model (SM). At 
its design center-of-mass energy, new physics production 
cross sections may be significant so that a data sample 
of modest integrated luminosity, 100 pb~^ or less, may 
contain a large number of new particles. The challenge is 
to distinguish events with new particles from those, many 
orders of magnitude more copious, attributed to the SM 
with limited understanding of the SM production rates 
and detector performance in early LHC data. 

Missing transverse energy, I^t, has discriminating 
power to reveal new particles interacting weakly with or- 
dinary matter produced via high energy parton collisions 
in laboratory conditions at the LHC. These weakly in- 
teracting particles may comprise the dark matter of our 
Universe [l|. They are expected in new physics mod- 
els, such as Pt^jarity conserving super-symmetry ^ and 
many others [3| . Missing transverse energy allows to per- 
form a broad search sensitive to the presence of such 
particles in collision data and is an observable that may 
lead to an early discovery at the LHC At the same 
time, missing transverse energy is one of the most dif- 
ficult observables to measure precisely and simulate ac- 
curately @ because it is measured by multiple detector 
sub-systems and subject to mis-measurements and back- 
grounds in any of them. 

In this paper, I discuss a new method 6] to predict 
backgrounds at high I^t for new physics searches in 
signatures consistent with SM V^-|-jets and ti+]ets 0, 
where F is a Z, 7 or MK. I assume that new particles 
are heavy and decay to SM particles emitting multiple 
jets so that high sensitivity is expected at high I^t and 
a large number of jets. Since main sources of artificial 
I^T come from the system of jets, the detector and non- 
collision effects, I model the instrumental response to the 
system of jets in t^-|-jets and other effects at high I^t in- 
situ using multi-jet QCD events. This method comple- 
ments and extends the work of Ref. where events 
with high rapidity objects arc used to model SM F-|-jets 
and multi-jet backgrounds in new physics searches with- 
out heavy reliance on IpT- The emphasis of this work, as 
that of Ref. Q , is on robustness against imperfections of 



background modeling required for new physics searches 
in early LHC data. 



II. OVERVIEW 

Monte Carlo (MC) simulation capable of modeling the 
detector response to SM processes is a great asset in new 
physics searches. However, there are two challenges in 
searches of early LHC data based on MC simulation. 
First, the SM F-|-jets production rates are difficult to 
predict from first principles. MC techniques are unre- 
liable in predicting backgrounds with a large number of 
jets and need to be tuned with high ^/s data. Theory cal- 
culations at sufficiently high order in many cases do not 
exist 9]. The structure functions have significant uncer- 
tainties in the small x range accessible at the LHC [lo| . 
Second, significant uncertainties in the calibration of the 
experimental apparatus are expected in early data tak- 
ing. Missing transverse energy is an observable that is 
particularly difficult to measure precisely and simulate 
accurately, since large jet energy fluctuations, detector 
artifacts, collision related and non-collision effects can 
produce non-Gaussian high I^t tails. These artificial 
l^T tails may resemble a signature of a new weakly in- 
teracting particle. 

To introduce the method, let us consider an event with 
a Z boson reconstructed in the di-muon channel and four 
jets. The four-momentum of the Z is well-measured so 
that the system of the four jets and other effects unre- 
lated to the di-muon system are the main source of I^t in 
this event. In order to develop a search in I/^t based on 
MC simulation, one would need to identify, understand 
and simulate the detector response to each of these ef- 
fects. Instead, I model these effects in-situ using multi- 
jet QCD events as follows. A sample of QCD events with 
four jets that have approximately the same configuration 
as the four jet system of the Z+jets event is selected. A 
l^T prediction, or a template, for this Z-|-jets event is ob- 
tained using the I^t distribution measured in the selected 
QCD sample and normalized to unity. This procedure is 
repeated for all other Z-|-jets events. The I^t templates 
are summed up to obtain a SM I^t prediction for the 
entire Z-|-jets sample. 

The photon momentum in 7-l-jets is also well- 
measured so that the same procedure applies to ob- 
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tain a SM I^t prediction in the 7+jets sample. In 
VF+jets and t<+jets with one of the two top quarks de- 
caying semileptonically, the ^+jets+^T signature, there 
is genuine Ifr from the undetected neutrino produced 
in W decays. To avoid reliance on MC and theory, I 
model the neutrino \pt\ spectra using charged lepton 
\'Pt I spectra. If W bosons are not polarized in the trans- 
verse plane, the two spectra should be the same. Event 
selection and polarized W bosons produced in top quark 
decays lead to differences in the charged lepton and neu- 
trino spectra. However, these differences are small and 
can be accounted for by corrections. A prediction for ar- 
tificial IpT hi W^+jets and tt-\-]ets is derived in the same 
manner as for Z+jets and combined with a modeled neu- 
trino \pt I spectrum to obtain a SM I^t prediction in the 
Z-f jets+^T final state. 

Since higher sensitivity to new physics is expected in 
events with a large number of jets, the focus is to model 
high I^T region in l^-l-jets events with 3 or more jets. 
Events with 2 jets are valuable as a validation and cal- 
ibration sample. This method is developed for searches 
in early LHC data. It will work best if the LHC start- 
up is quick, new particles are strongly produced and not 
very heavy, e.g., such as squarks and gluinos in the low 
mass mSUGRA CMS and ATLAS benchmarks [6]. With 
this in mind, a prediction of SM backgrounds in high 
Ifx tails to about 20% may be sufficient to reveal new 
physics. For this reason, an accuracy benchmark for this 
method is to predict SM backgrounds in V-l-jets at high 
Ifx and a large number of jets (3 jets or more) to about 
20% or better in a data-driven manner. 



III. EXPERIMENTAL ASPECTS 

The CMS and ATLAS experiments use multi-purpose 
detectors at the European Organization for Nuclear Re- 
search (CERN). Detailed descriptions of the detectors 
can be found in Ref. The detectors are capable of 

reconstructing electrons and muons with high efficien- 
cies and low fake rates for lepton \pt\ > 20 GeV in the 
\r]\ < 2.5 range (In this paper the symbol I is used 
to denote an e, /i, but not r. Charge-conjugate modes 
are implied.) In both detectors, photons and jets can 
be reconstructed reliably within \ri\ < 2.5 and \ri\ < 3.0, 
respectively. 

To study the method, mock data samples are gener- 
ated for the following SM processes: Z-l-jets (5.0 fb~^, 
up to 5 partons, Z — + W^H-jets (1.0 fb~-^, up to 

5 partons, W — > lui), tt-fjets (1.0 fb~^, up to 4 par- 
tons, tt — > Ivibbjj), 7-fjets (400.0 pb^^, up to 5 
tons) and QCD jets (1.0 pb~^, up to 5 partons) 
(The same samples are used in Ref. Q.) The inte- 
grated luminosity listed in parentheses is used every- 
where in tests in this paper but section IxH These sam- 
ples are generated with ALPGEN [3] at the parton level. 
PYTHIA [3] is used for parton showering, hadroniza- 
tion, simulation of the underlying event and jet recon- 



struction. To model features of a new physics signal in 
search distributions, mock signal data samples for Min- 
imal Supergravity (mSUGRA) benchmark points LMl 
and LM4 [3,[il,[l3 are generated with PYTHIA. 

Electrons and muons are required to have \pt\ of at 
least 20 GcV in the \ri\ < 2.5 range. Photons are re- 
constructed above the \pt\ threshold of 30 GeV in the 
\ri\ < 2.5 range. Jets are reconstructed using the PY- 
CELL algorithm [l^ and required to be within \r]\ < 3.0. 
A low jet \pt\ threshold of 20 GeV is used in the I^t mea- 
surements in order to collect the energy deposited in the 
calorimeters to a fuller extent. Higher jet \pt\ thresholds, 
50 GeV or more, are used to measure other observables 
in a robust manner as indicated below. I require that the 
leading jet and I^t be not aligned in the transverse plane 
within 0.15 radians: 0.15 < | A^i<=ad jot-^T| < (7^ _ 0.15). 
(The jet with the highest \pt\ in an event is the leading 
jet of this event. Any other jet in this event is a non- 
leading jet.) It is assumed that the triggering and event 
reconstruction efficiency in each channel is 50%. 



"0.45 

^-0.4 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 






— Total resolution 


\ \ 
\ \ 


Stochastic term 

— Noise term 

- - Saturation term 







10^ 



10" 
[GeV] 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



pReco'd ^ pTrue 



FIG. 1: Top: the Gaussian contribution to the jet energy 
resolution as a function of true jet jprl is shown in the solid 
line. The noise, stochastic and saturation contributions to 
the jet resolution function are shown separately. Bottom: the 
jet energy smearing functions for 500, 250, 100 and 50 GeV 
\Pt\ jets are shown in the solid, dashed, dot-dashed and dot- 
ted lines. 

The vector of missing transverse energy, I^t, is calcu- 
lated as the vector opposite to the sum of pr measure- 
ments of charged leptons, photons, and jets in each event. 
Since electron, muon and photon momenta are measured 
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with high precision compared to jets, their contribution 
to artificial Ifx in events with a large number of jets is 
negligible. To emulate detector resolution effects for jets, 
jet energies measured by PYTHIA are smeared. The jet 
smearing function has three components: a) a Gaussian 
with 

ct{\pt\)I\pt\ = V (7-0/IptI)^ + (1.2//^)2 + (0.04)2, 

where \pt\ is measured in GeV, b) an exponential low- 
side tail that stretches from [1.0 — 2cr(|pT|)/|PT|] to 0.0 
added to the Gaussian component 3% of the time, and 
c) similarly, an exponential high-side tail from [1.0 + 
2o'(|pt|)/|pt|] to infinity added to the Gaussian 1% of 
the time. Figure [T] shows the \pt \ dependence of the 
Gaussian smearing and the full smearing function with 
the non-Gaussian tails for a few fixed jet \pt\ values. 
This jet smearing function is constructed based on stud- 
ies of the CMS and ATLAS detectors [1, ^ to represent 
the jet response characteristic of the two detectors. Since 
the jet system tends to be the dominant source of artifi- 
cial all three components of the jet smearing function 
are varied in the studies of robustness and limitations of 
the method discussed below. 

The selection criteria used in the paper are not opti- 
mized to any new physics model. Instead, they are cho- 
sen to ensure robust detector performance and maintain 
sensitivity to a wide range of new physics models at high 
Ij^T and a large number of jets. The selection criteria 
can be modified without producing a significant effect on 
the method's performance. The mSUGRA benchmarks 
listed above are used only to illustrate generic features 
of a new physics contribution. The goal of this paper is 
to demonstrate the scope of the method and its perfor- 
mance rather than to attain high sensitivity to a specific 
model and quantify it. 



IV. ALGORITHM 

I will describe in detail the algorithm in this section 
and present results of its tests in the next section. Let 
us assume that the V momentum in a y-|-jets event is 
known. The resolution and other effects producing artifi- 
cial missing transverse energy for this event are modeled 
using multi-jet QCD events with a kinematic configura- 
tion of jets similar to that in the F-|-jets event. 

A prominent difference between the jet systems in 
F-l-jets and QCD is that in F-|-jets events the jet sys- 
tem recoils against the V while in QCD events it is at 
rest in the transverse plane. A key is to select multi-jet 
QCD events to predict I/^t in a manner that captures 
effects generating artificial Ifr in-situ but allows for the 
difference stemming from the boost of the jet system in 
the transverse plane in y-|-jets. 

It is not necessary to model accurately every degree of 
freedom in the jet system of l^-|-jets on an event-by-event 
basis by QCD for two reasons. First, each l/-|-jets event 



is modeled by a large sample of QCD events so that mis- 
modeled correlations are averaged out over this sample 
of QCD events. Second, I^t is measured for the entire 
V^-|-jets sample so that mis-modeled correlations are av- 
eraged out over V^-|-jets events as well. These averaging 
effects allow to develop a simple algorithm. 

Multi-jet QCD events are selected using two variables: 
(1) Nj, number of jets above a high \pt\ threshold 
(50 GeV or higher), and (2) Jt = \Pt^\ fo^' j^ts above 
a low 20 GeV \pt\ threshold (the same jet threshold is 
used for I^t measurements) [l3|. A QCD sample is se- 
lected for each pair {Nj, Jt bin), the width of Jt bins is 
10 GeV (100 GeV) below (above) 1 TeV. A template 
is obtained for each of these samples as a Ifr distribu- 
tion in that sample normalized to unity. For each y-|-jets 
event, Nj and Jt are measured and used to select the 
I^T template with the same Nj in the corresponding Jt 
bin, which represents an artificial I^t prediction for this 
V^-fjets event. These templates are summed up over all 
V^+jets events to obtain a I^t prediction for the entire 
V+]ets sample. 

Two sets of jet \pt\ thresholds are used to measure Nj. 
In the first set, the jet threshold for Nj is 50 GeV. 
QCD events for this jet threshold can be collected us- 
ing prescaled low jet triggers. In the second set, 
the jet thresholds are equal to (or higher than) the jet 
\pt\ thresholds that can be used in unprescaled multi-jet 
triggers. For the second set, I use 140 GeV for Nj — 2, 
80 GeV for A^j = 3 and 60 GeV for Nj > 4. These jet 
\pt I thresholds can be changed depending on the trigger 
rates in data without significant effect. 

The V momentum in Z-|-jets and 7-l-jets is well- 
measured so that the application of the algorithm is 
straightforward in these channels. A comparison of pre- 
dicted and observed yields is shown in Figures [5] and [3] for 
Z-|-jets and 7-l-jets, respectively, where the top (bottom) 
row shows results for the 50 GeV (high) jet \pt\ thresh- 
olds for Nj. It is seen that the prediction (solid line) 
is very good for Nj = 3 and Nj > 4. For events with 
Nj = 2, the measurement (dashed line) is about 20% 
below the prediction at I^t < 20 GeV. The mechanism 
responsible for this bias is discussed in section |Vll Since 
new physics is not expected to contribute at small I^t, 
to remove this bias, the prediction is normalized to the 
measurement in the I^t& [50, 100] GeV interval. This is 
done for the I^t predictions in Z-|-jets and 7-t-jets events 
with Nj = 2 everywhere in the rest of the paper. 

In M^-|-jets and tt-fjets events reconstructed in the 
Z-|-jets-|--^T channel, there is genuine missing transverse 
energy from undetected neutrinos produced in W decays. 
To study only the effect of ]^t mis-measurements, I con- 
sider the dominant W —^ Ivi and tt — > Ivibbjj contribu- 
tions and assume that the neutrino spectra are known 
until section IIXI To model I^t resolution effects, the 
neutrino \pt \ 's are smeared with the artificial I^t predic- 
tions obtained from multi-jet QCD events. This is done 
on an event-by-event basis assuming that the neutrino 
Pt and the artificial I^t interfere at at a random an- 
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FIG. 2: Algorithm performance in Z+jets for Nj = 2, 3 and > 4 in the first, middle and third columns, respectively. The top 
row shows results for the 50 GeV jet \pt\ threshold for Nj. The bottom row is for the high jet \pt\ thresholds for Nj. The 
observed distributions are shown in the dashed lines, their predictions obtained using multi-jet QCD events are the solid 
lines. 
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FIG. 3: Algorithm performance in 7+jets for Nj = 2, 3 and > 4 in the first, middle and third columns. The top (bottom) 
row shows results for the 50 GeV (high) jet Iprl threshold(s) for Nj. The observed distributions are shown in the dashed 
lines, their predictions are the solid lines. 
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FIG. 4: Algorithm performance in Z-|-jets-|-_^T for Nj = 2, 3 and > 4 in the first, middle and third columns. The top (bottom) 
row shows results for the 50 GeV (high) jet \pt\ threshold(s) for Nj. The observed distributions are shown in the dashed 
lines, their predictions are the solid lines. 



gle (j) distributed uniformly from to tt in the transverse plane. Figure [4] shows how well the method works for the 
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50 GeV (top row) and high (bottom row) jet thresh- 
olds for Nj in the l+jets+l^T final state, where both 
VF+jets and tt+jets are included according to their ex- 
pected production cross sections. 

In section IIXI it is demonstrated that one can ap- 
proximate the neutrino \pt\ spectra by the charged lep- 
ton \pt\ spectra. The contribution from W — > tVt in 
M^-l-jets and tt-|-jets is also considered in section IIXI 
With these extensions, the method can be used to predict 
the l^T distribution in the ^-|-jets-|-.^T final state, which 
has high sensitivity to a variety of new physics models 
with new weakly interacting particles in early data. 
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FIG. 5: Ratios of observed and estimated integrated yields for 
Z-j-jets (top), 7-|-jets (middle), Z-j-jets-j-^r (bottom) obtained 
for the 50 GeV jet \pt\ threshold for Nj. In each plot three 
types of markers are shown for Nj — 2 (circles), 3 (squares) 
and > 4 (triangles). The shaded markers for Z-j-jets and 
7-|-jets show the ratios before the predictions are normalized 
at low as described in the text. Note, the ratios are 
correlated since yields are integrated upwards. 

For brevity, in the rest of the paper, I present results 
of studies for the 50 GeV jet threshold used to measure 
Nj. They have higher statistical precision than those for 
higher jet \pt\ thresholds for Nj. Ratios of observed 
and predicted yields, A^obscrvcd/^Piodictcd, are shown in 



Figure [5l where the yields in each I^t bin are integrals 
of the distributions shown in the top rows of Figures [2] 
through d] from that bin's I^t value to infinity. The 
algorithm performs at least as well when the set of higher 
jet \pt\ thresholds for Nj is used. 

Since the QCD production cross section is very large 
at the LHC, only a small QCD sample is needed for this 
method to work, e.g., I pb^^ of QCD is used to model 
I^T distributions in 5 fb^-'^ of Z-|-jets in this paper. 
Again, the QCD sample for templates can be collected 
via prescaled small \pt\ jet triggers and unprescaled high 
\pt\ multi-jet triggers. Due to the large QCD production 
cross section, the relative contribution from electroweak 
processes with genuine I^t from neutrinos in this sample 
is negligible for searches in early data. 



V. ROBUSTNESS 

The goal of this method is to capture effects generat- 
ing high artificial I^t in-situ using multi-jet QCD events. 
To demonstrate how well the method works, I present a 
set of tests in which increased jet mis-reconstruction is 
introduced. In each test, an identical change to the mock 
data samples for F-l-jets and the QCD sample is made 
and the analysis procedure is repeated. Figure [6] shows 
how drastic the effect of these changes on the I^t dis- 
tribution in 7-|-jets can be. (The 7-l-jets channel in 
Figure [6] is used since there is no genuine I^t in this final 
state and it has a larger yield than Z-|-jets.) Test results 
are presented in Figure [7] for Nj = 3 as ratios of observed 
and estimated integrated yields. For brevity, test results 
for Nj = 2 and > 4 are not shown but discussed and 
compared to those in Figure [T] 

Jet reconstruction efficiencies are not equal to unity. 
To test if the method models effects due to undetected 
jets accurately, an identical source of jet inefficiency is 
introduced in y-|-jets and QCD events. I remove jets 
that faU in a veto cone of AR = Arf + A^^ ^ 0.8 at 
{•q,4>) — (0.0,0.0) where is the azimuthal angle. 
Since softer jets are more likely to be lost, only non- 
leading jets are removed in the veto cone in this test. 
The effect of this inefficiency on the I^t distribution in 
7-fjets for Nj = 3 can be assessed by comparing the 
solid line with the dashed grey line in Figure [6l Test 
results for Z-|-jets, 7-l-jets and /-|-jets-|--^T, W^+jets and 
tt-|-jets, with Nj = 3 are shown in Figure [7] in circular 
markers. The increased artificial I^t tail due to lost jets is 
modeled accurately by the method. The same conclusion 
holds for iVj = 2 and > 4. 

In the next test, the jet energy smearing is increased. 
Two tests are made: (a) the Gaussian cr (|pV|) is dou- 
bled with the area and shape of the non-Gaussian tails 
unchanged, and (b) the area of the low and high-side 
non-Gaussian tails is doubled with the Gaussian compo- 
nent kept unchanged. The effect of the additional jet 
energy smearing in test (a) on the I^t distribution in 
7-l-jets for Nj = 3 is shows in Figure [S] in the dotted 
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FIG. 6; Illustration of effects associated with jet mis- 
reconstruction on artificial in 7+jets in the Nj — 3 
bin for the jet \pt\ threshold of 50 GeV. The black line is 
a reference distribution from Figure [3] Jet reconstruc- 
tion inefficiencies (dashed grey), increased jet energy smear- 
ing (dotted grey) and extraneous energy (dot-dashed black 
and dotted black) from the tests in section |V] and |Vl] signifi- 
cantly increase artificial ^t- 



grey line. Ratios of observed and estimated yields from 
test (a) are in square markers in Figure [7] for Nj = 3. 
Similar results are observed in the other two Nj bins, and 
in test (b). Again, templates constructed from multi-jet 
QCD events capture effects from additional jet smearing 
in-situ so that the level of background at high I^t in 
F-fjets is predicted accurately. 

Hot cells or noise in the calorimeters, backgrounds 
from the proton beams, cosmic rays, underlying event 
or additional pp interactions in the same bunch crossing 
contribute extra energy and jets erroneously attributed 
to those produced in T^-|-jets and QCD processes. Since 
additional jets have higher probability to be soft, I test 
the method's ability to model such effects by adding ex- 
tra jets with a soft uniform spectrum from to 
50 GeV with a 20% probabihty to each T^-|-jets or QCD 
event. These extra jets change Jt and I^t, but do not 
change Nj. I find that all predictions are good in this 
test. 

I repeat the previous test with a uniform \pT \ spectrum 
of additional energy contributions covering the range 
from to 100 GeV added with a 10% probability to 
y+jets and QCD events. This produces a strong effect 
on the I^T distribution shown for 7-l-jets with Nj = 3 in 
the dot-dashed line in Figure (S] Ratios of observed and 
estimated yields for the three V-l-jets processes in the 
Nj = 3 bin are in triangular (down) markers in Figure [71 
I find that the prediction is robust in the Nj = 3 and 
> 4 bins. In the Nj = 2 bin in Z-|-jets and 7-l-jets, a 
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FIG. 7: Ratios of observed and estimated integrated yields 
in Z-l-jets (top), 7-l-jets (middle), VK-fjets and tt-|-jets (bot- 
tom) all for Nj = 3 and the 50 GeV jet \pt\ threshold for Nj 
from robustness tests in section|V] Circles, squares, triangles- 
down, triangles-up are for tests with increased inefficiencies 
for non-leading jets, increased Gaussian jet energy smearing, 
extraneous energy contributions and a modified Nj compo- 
sition of the QCD sample. Note, ratios in the plots are cor- 
related as tests are made using events drawn from the same 
mock data samples, and yields are integrated upwards. 



bias is observed. The origin of this bias stems from dif- 
ferences in Nj and Jt distributions between F-|-jets and 
QCD. I discuss it and more stringent tests with extrane- 
ous energy contributions in the next section. 

The cross section ratios for F-l-jets and QCD pro- 
cesses: a^+^^^^n jets)/ cr<^™(n jets), (7^+j*='"(n jets)/ 

^V+}cts^^ -I- 1 jets) or crQCD(-^ jg^g^ crQCD(^ _^ jg^g)^ 

where n is equal to 2 or more, in LHC data are likely to 
differ from that of ALPGEN used in this study. There 
may also be differences in other differential distributions 
in the jet system of y-fjets or multi-jet QCD events be- 
tween LHC data and ALPGEN. To test how sensitive 
the method is to such differences, I vary the ALPGEN 
ratios a'^'^^{n jets) / a'^'^^ {n + 1 jets), n > 2, by a fac- 
tor of 1.5 up or down. Test results with reduced ratios 
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for Nj = 3 are shown in Figure [7] in triangular (up) 
markers. The I^t predictions are good in this test be- 
cause they are made on an event-by-event basis using 
QCD events with the same Nj and Jt- QCD events 
with other values of Nj and Jt are included only if they 
are misreconstructed, which is a second order effect, but 
it can become significant in regimes where distributions 
fall or rise steeply. Test results for Nj = 3 and > 4 are 
all good. For Nj = 2 in Z-|-jets and 7-l-jets, when the 
cr^^^(n jets)/cr^^^(n -|- 1 jets) ratios are reduced the 
prediction improves; when the ratios are increased the 
prediction becomes biased. The origin of this bias is the 
same as in the previous test and is discussed in the next 
section. 

In conclusion of this section, the quality of the I^t pre- 
diction improves at larger Nj. The I^t prediction is 
robust for Nj — 3 and > 4 in all tests. Events with 
Nj = 2 are more susceptible to biases for two reasons. 
First, there are significant differences in the differential 
distributions describing jets in QCD and F-|-jets: in QCD 
di-jets, the jets come mainly from leading order parton 
interactions, while in V^-|-jets, the jets are from higher or- 
der processes. Second, the averaging effects discussed in 
section IIVI are not as strong when the number of jets is 
small. Nevertheless, only two tests for Nj — 2 are biased 
in this section. Any other effect that generates artificial 
I^T in the same manner in the jet system of y-|-jets and 
multi-jet QCD events should be modeled in-situ by the 
method. I next discuss the method's limitations revealed 
in more stringent tests. 



VI. LIMITATIONS 

I increase the degree of jet mis-reconstruction up to a 
point where the method becomes biased to explore the 
boundaries of the domain where the method works. This 
allows to understand in greater detail mechanisms that 
may lead to a bias. At the end of this section, I discuss 
how to avoid regimes where the method is biased. 

The test with the jet veto cone introduced in the pre- 
vious section is repeated with a modification such that 
leading jets falling into the veto cone are removed. This 
is a stringent test since leading jets are less likely to be 
undetected. Test results are shown in Figure [S] in circu- 
lar markers, for brevity, only for 7-l-jets in the Nj — 3 
bin. While the prediction partly takes into account the 
effect of undetected leading jets, it underestimates the 
background at high I^t in that Nj bin in Z-|-jets and 
7-l-jets. The prediction is biased because in QCD events 
I^T is always less than Jt, by the definition of I^t in 
section IIIII In l^-l-jets, y is a Z or 7 here, I^t can be 
greater than Jt when the leading jet recoiling against an 
energetic V boson in the transverse plane is lost. The 
y-|-jets events with I^t larger than Jt can not be mod- 
eled by the algorithm in section llVl This bias is larger for 
Nj = 2, while in the Nj > 4 bin, the prediction is good 
for both Z-l-jets and 7-l-jets. In l+^ets+l^T, W+iets and 



tt-|-jets combined, due to a genuine I^t contribution from 
neutrinos to the full I^t, this bias does not appear. 
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FIG. 8; Ratios of observed and estimated integrated yields in 
7-l-jets for Nj = 3 and the 50 GeV jet \pt\ threshold for Nj 
from tests in section IVII Circles, squares, triangles-up and 
triangles-down are for tests with increased inefficiencies for 
leading jets, increased non-Gaussian jet energy fluctuations 
and extraneous energy contributing in jets or anywhere in 
the calorimeters. Note, these ratios are correlated as tests are 
made using events drawn from the same mock data samples, 
and yields are integrated upwards. 

I repeat the test with increased jet energy mis- 
measurements after tripling the area of the lower non- 
Gaussian tail in the jet response function and reducing 
the magnitude of its slope on the logarithmic scale of the 
lower plot in Figure [T] by a factor of 2. The prediction is 
biased in the Nj = 2 bin for both Z-|-jets and 7-l-jets. 
The quality of the prediction improves in the Nj = 3 
bin, shown in Figure [8] in square markers for 7-l-jets, and 
it is good for Nj > 4 in Z-|-jets and 7-l-jets. One should 
expect a bias for large low-side tails in the jet response 
function appearing via the same mechanism operating 
in the previous test. The ^-l-jets events containing jets 
fiuctuated down in \pt\ can have Jt that is less than 
I^t- Such events can not be modeled by the algorithm of 
section ITVl In the Z-|-jets+^T final state, the prediction 
of the full I^t, which includes the neutrino momentum, 
is good in all Nj bins. Note, for a large low-side tail in 
the jet response function, the jet energy scale may be- 
come biased. Effects due to a jet energy scale offset are 
discussed in section IVIIII 

I make two tests with a harder spectrum of addi- 
tional energy contributions unrelated to 7+jets and 
QCD events. In these tests, the spectrum of additional 
energy is uniform in Et from to 1 TeV added with a 
1% probability to both V^-|-jets and QCD. Since there is 
no genuine I^t in 7-l-jets and QCD, the requirement on 
|^0ioadjet-|;T| (section Uni that and the leading jet 
be not aligned in the transverse plane removes a fraction 
of events with high Et extraneous contributions. In the 
first test, additional energy depositions contribute only 
to jets that are above the \pt\ threshold for Nj. The 
I^T distribution in 7-l-jets, Nj = 3, is shown in the dotted 
black line in Figure [S] with a large artificial high I^t tail. 
Ratios of observed and estimated yields are in triangu- 
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lar (up) markers in Figure [8] In 7+jets, the prediction is 
good for Nj — 3 and > 4, and it is biased in the Nj — 2 
bin for the following reason. The Jt spectrum in QCD 
events tends to be softer than that in V+jets events 
with the same Nj. (This effect is most pronounced for 
Nj = 2.) The fraction of soft QCD multi-jet events pro- 
moted to higher Jt by extraneous energy depositions 
tends to be larger than that fraction in IZ-l-jets. Since 
such events have larger I^t due to the extraneous energy 
depositions unbalanced in the transverse plane, the level 
of background at high I^t is overestimated. 

In the second test, extraneous energy contributions 
are added randomly in the rj — (jj plane so that Nj also 
tends to increase. Ratios of observed and estimated 
yields for 7-t-jets, Nj — 3, are in triangular (down) 
markers in Figure [51 The prediction overestimates the 
background in all Nj bins. This happens because 
ctqcd(« jets)/(TQCD(n + 1 jets), n > 2, is higher than 
cry+jets("- jets) / (Jv+]cts{'n+ 1 jcts) in the mock data sam- 
ples so that the fraction of events with Nj — n recon- 
structed erroneously in the Nj — {n + 1) bin due to 
an extra energy deposition is higher in QCD compared 
to that fraction in V"+jets. (Again, this effect is most 
pronounced for Nj — 2.) Since these mis-reconstructed 
events have larger I^t, the prediction overestimates the 
background. 

Test results with extraneous energy contributions for 
Z+}ets are qualitatively similar to those for 7-l-jets. In 
?-|-jets-|-^T, the biasing effects discussed above are inter- 
twined with additional effects due to the presence of a 
neutrino in the final state and the tt-l-jets contribution. 
The genuine I^t from the neutrino makes the require- 
ment on I Af/)''^'*'^ less efficient in suppressing high 
Et extraneous contributions. The tt-|-jets events con- 
tribute to further differences in Nj and Jt spectra be- 
tween W^-|-jets and QCD. I find that in /-|-jets-|--^T the 
prediction tends to overestimate the background in the 
tests with extraneous energy depositions, and the quality 
of the prediction improves with Nj. 

The explanation for biases with extraneous energy con- 
tributions can be generalized and applied to other sources 
of biases in the predictions of artificial I^t- Biases in 
this method can appear when two conditions are satis- 
fied. First, an observable 0^°^^ describing the jet sys- 
tem and related to I^t is severely mis-measured lead- 
ing to increased artificial I^t and mismeasured Nj or 
Jt- Second, there are significant differences in the dif- 
ferential distribution of O^^^^ between multi-jet QCD 
and F-|-jets events on the scale associated with the mis- 
reconstruction. When these two conditions are satisfied, 
the fraction of events with mis-measured Oj*^**^ is different 
in QCD and F-l-jets for the same reconstructed Nj and 
Jt- If the mis-measured fraction of multi-jet QCD events 
is larger, the prediction overestimates the background. 
Otherwise, the predicted background is underestimated. 

Regimes with severely misreconstructed events where 
the method may become biased need to be avoided. By 
imposing event quality criteria or improving the jet re- 



construction, e.g., using the tracking systems, [2^ one 
can reduce the number of such events. Moderately mis- 
reconstructed events are modeled in-situ by the method. 
The F-|-jets sample with Nj = 2 is the most challeng- 
ing for this method. The method performs better at 
higher Jt and Nj, where the sensitivity to new physics 
is higher. There are several reasons for that: a) there are 
fewer differences between the hadronic systems in y-|-jets 
and QCD, b) the averaging effects over F-|-jets and QCD 
events are stronger and c) the jet reconstruction performs 
better at higher jet |pt|- 



VII. tt+JETS 

SM tt-|-jets events, where tt — > Ivibbqq, constitute a 
dominant background in the l+}ets+l^T signature for 
Nj > 3. The shapes of Nj and Jt spectra in these 
events differ from those in QCD events collected for tem- 
plates and from those in F-l-jets. The calorimeter re- 
sponse to &-jets in tt-|-jets differs from that of light quark 
and gluon jets [2l|. These effects lead to a bias in the 
prediction of artificial I^t in tt-l-jets. To demonstrate 
this bias clearly. Figure [9] (top) shows the artificial I^t in 
t<+jets for Nj > 4 (dashed line), where the neutrino 
four-momentum is assumed to be measured so that it is 
included in the I^t calculation, and its prediction using 
QCD templates (solid line). Note, at large I^t, this bias 
is an order of magnitude smaller compared to the gen- 
uine I^T from the neutrino in the final state having the 
\pt\ spectrum shown in the dot-dashed line in the same 
Figure. When the neutrino \pt\ spectrum is combined 
with the artificial I^t in the full I^t prediction, the bias 
becomes insignificant as seen in Figure [H] (bottom) . 

The artificial I^t is a dominant contributor in events 
with small genuine I^t Figure [H] (bottom) shows that the 
accuracy of its prediction is sufficient to model the full 
I^T distribution at small I^t At high I^t, the missing 
momentum from the neutrino dominates over artificial 
I^T SO that the accuracy of the full I^t prediction is 
highly dependent on how well the neutrino spectrum is 
modeled. The modeling of neutrino spectra is discussed 
in section HXl 

Despite the fact that the bias in the artificial I^t pre- 
diction for tt-t-jets is insignificant in the full I^t pre- 
diction in Z-|-jets-|--^T, it is instructive to examine how it 
behaves when selection criteria or the algorithm of sec- 
tion IIVI are modified. Two observations can be made. 
First, the bias becomes smaller when the jet \pt\ thresh- 
old for I^T and Jt is reduced or the rj coverage for jets is 
increased since the total energy is collected to a fuller ex- 
tent with more inclusive requirements. Optimal require- 
ments on these variables can only be determined using 
data because at smaller \pt\ and larger \r]\ more noise 
and backgrounds are expected. Second, in tt-|-jets, there 
tends to be more jets included in the I^t and Jt calcula- 
tions that are below the jet |pV| threshold for Nj. Since 
the jet resolution improves as the jet \pt\ grows, the 
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FIG. 9: Top: artificial in tt+jets (dashed) and its pre- 
diction (solid) for Nj > 4 and the 50 GeV jet \pt\ threshold 
for Nj. Bottom: full (dashed) and its prediction (solid) 
combining both the neutrino \pt\ spectrum and the artificial 
prediction for the same Nj and jet \pt\ selection. In both 
plots the neutrino spectrum is shown as the dot-dashed line. 



prediction can be improved by making I^t templates in 
coarse bins of R{Jt) = J^^ /Jt, where J^'® is a scalar 
sum of jet \pt\'s for jets above the \pt\ threshold for 
Nj. Alternatively, the same effect can be achieved by 
modifying the composition of the QCD sample used for 
templates. 



VIII. JET ENERGY SCALE 

Jet energy measurements can be systematically biased 
in early data. Let us consider a case when jet ener- 
gies are under- measured uniformly in jet \pt\- Such mis- 
measurements cancel to first order in I^t measurements 
in QCD events, fn F-|-jets, since the jet system recoils 
against the V, the jet energy mis-measurements add up 
coherently along the V direction in the transverse plane. 
To avoid a bias due to this difference, the jet energy scale 
needs to be calibrated. Since the method is capable to 
model large tails in the jet response function, a precise 
calibration of the jet energy scale as a function of rj and 
(j) (azimuthal angle) is not required. The jet energy scale 
can be calibrated with sufficient accuracy using standard 
techniques based on 7-l-jets and Z-j-jets {Nj < 2) pro- 
cesses [20] in very early data. 

Figure (TU] gives a comparison of I^t distributions in 
Z-|-jets (top) and Z-|-jets-|--^T (bottom) for Nj — 3 with- 
out (solid) and with a 10% (dashed) and 20% (dotted) jet 



FIG. 10: The distributions without (solid) and with 
20% (dotted) or 10% (dashed) jet energy scale offset in 
^-fjets (top) and Z-|-jets-|-_^T (bottom) for Nj — 3 with the 
50 GeV jet threshold for Nj. The distributions in each plot 
are normalized to the same area. 



energy scale offset downwards uniform in jet \pt\- The 
prediction becomes good for a 10% or smaller offset in 
Z+jets and 7-l-jets. One may reduce the effect from a 
residual jet energy scale offset on I^t in Z+jets and 
7+jets by normalizing the predicted I^t shape to the 
observed distribution in the small I^t region, for exam- 
ple, for l^T G [50, 100] GeV. Demands on the precision 
of the jet energy calibration in Z+jets+^r are higher. 
Finally, even before the jet energy scale is calibrated, 
one can make a search in the projection of I^t on the 
axis perpendicular to the V direction (the I direction in 
Z+jets+^r) in the transverse plane, I^tt- Searches in 
l^TT are less sensitive to effects associated with the jet 
energy scale offset since those lead to a bias along the 
V direction. 



IX. NEUTRINO SPECTRA IN W DECAYS 

In the I + jets+^T signature, dominated by 
T^+jets and tt+jets, there are one or more undetected 
neutrinos in the final state. To model I^t in these events, 
one needs a prediction or a measurement of the neutrino 
\pt\ spectra, which can be combined with I^t reso- 
lution predictions from QCD templates. The neutrino 
\pt\ spectra could be obtained from MC simulation. Or, 
the neutrino \pt\ spectra can be modeled in a data- 
driven manner using charged lepton ]pV] spectra as de- 
scribed in this section. 
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The solid and dashed Unas in plot (a) of Figure [TT] are 
the neutrino and charged lepton spectra in VF+jets for 
Nj = 3, W ^ li'i, where the requirement on the charged 
lepton \pt\ of at least 20 GeV is removed. Since 
W bosons are not polarized in the transverse plane in 
W^+jets, the neutrino and charged lepton \pt\ spectra 
are consistent. The application of a \pt\ threshold on 
the charged lepton makes its spectrum harder, while the 
neutrino spectrum becomes softer, as seen in plot (b) of 
the same Figure. To model the neutrino spectrum using 
the reconstructed charged lepton spectrum in H^+jets, 
the effect of the charged lepton \pt\ threshold needs to 
be corrected for. The corrections need to be obtained 
from MC simulation. 

The solid and dashed lines in plot (c) of Figure [TT] are 
the neutrino and charged lepton spectra in tt+jets for 
Nj = 3, — > Ivibbqq, without a threshold requirement 
on the charged lepton \pt\- In the SM, 30% of bosons 
in top decays are produced in the transverse-minus he- 
licity state (left-handed) and the rest are longitudinally 
polarized [23|. Left-handed bosons tend to produce 
charged leptons with a \pT \ spectrum that is softer com- 
pared to the neutrino \pt\ spectrum as seen in plot (c). 
Since the two spectra have similar shapes, it is possible to 
use the muon spectrum to model the neutrino spectrum 
in tt+jets. Again, when a charged lepton \pt\ threshold 
is applied, the charged lepton spectrum becomes harder 
while the neutrino spectrum becomes softer, which leads 
to a higher consistency between the two spectra seen in 
plot (d). Nevertheless, the effects of the W polarization 
in top decays and the event selection, mainly due to the 
charged lepton threshold, in tt-|-jets, in general, need 
to be corrected. 

In order to determine corrections to the charged lepton 
spectra for VF-|-jets and tt+iets from MC simulation, 
one needs to measure the shape of the |pV | dependence 
of lepton reconstruction efficiencies and the relative frac- 
tions of Vl^-l-jets and tt-fjets in the data sample. The for- 
mer can be readily done via a standard technique based 
on Z — > decays [l^. The latter should come from an 
independent measurement. With these two ingredients, 
corrections can be determined from MC simulation. 

Since corrections to the charged lepton spectra are 
small, the reliance on details of MC simulation to de- 
termine the neutrino \pt\ spectra is minimized. For a 
20 or 15 GeV threshold on charged lepton \pt\, no cor- 
rections are required to predict the I^t distributions in 
tt-|-jets in all Nj bins to 20% or better in the mock data 
samples. Corrections are needed for W-t-jets. The lower 
plot in Figure [TT] shows the I^t distribution and its pre- 
diction in /+jcts+-^T, W^-l-jets and tt-|-jets combined, for 
Nj > 4 based on the charged lepton spectrum without 
corrections. Since tt-|-jets dominates over W-l-jets in 
the Nj > 4 {Nj = 3) bin, the prediction is good to 
15% (25%) at high without corrections. The Nj > i 
bin, where the prediction is the most robust, is likely to 
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FIG. 11: Comparisons of charged lepton and neutrino spectra 
in VK-|-jets (a,b) and tt+jets (c,d) without (a,c) and with (b,d) 
charged lepton \pt\ thresholds all for Nj — 3. The lower plot 
shows a prediction based on the charged lepton spectrum 
in Z+jets-|-_^T, Nj > 4, without any corrections applied. In 
all plots, the threshold for Nj is 50 GeV. 



have the highest sensitivity to a new physics contribution 
compared to lower jet multiplicity events. 



B. 



W 



In the /-|-jets-f signature, there is background from 
tauonic W decays in W-|-jets and tt-|-jets. Tauonic W 
decays produce at least one additional neutrino that is 
a source of differences between the muon and neutrino 
spectra. 

There are two types of tauonic W decays that con- 
tribute significant background: (1) T^/'-l-jets and ti-|-jets. 
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where W~ TVr with r IviVr, and (2) ti events, 
where W~ — > Ivi and fv^. with t — > (hadrons Vr)- 

The contribution from tauonic W decays is an order of 
magnitude smaller compared to that from W —> Ivi de- 
cays. (The tauonic background of type 2 can be sup- 
pressed by vetoing events with isolated single hadronic 
tracks.) The r branching fractions are well known. 
Therefore, the effects from W TVr on IpT predictions 
can be well-modeled by an additional smooth correction 
to the muon jp^ | spectra that can be determined by MC 
simulation. Since the contribution from tauonic W de- 
cays is smaller compared to that from W Ivi decays in 
Z-|-jets-|-^T, corrections for them are not discussed further 
in this paper. 



X. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties need to account for the sta- 
tistical precision and biases of the method's background 
predictions at high I^t- Mechanisms by which biases may 
appear are discussed in section lVTl The method's suscep- 
tibility to them can be studied in both data and MC. 

A sample of events with Nj — 2 for a small jet 
\pt\ threshold, such as 50 GeV, is more sensitive to bi- 
ases. The relative contribution from new physics can not 
be large in this sample. Therefore, these events can be 
used to validate the method's performance and place an 
upper bound on its biases at higher Nj and jet \pt\- 
Similarly, the application of event quality criteria are ex- 
pected to reduce the number of severely misreconstructed 
events that may lead to biases in the prediction of arti- 
ficial I^T- By varying the event quality selection criteria, 
one can determine if the method is subject to such biases 
or estimate their size. An excess due to a new physics 
contribution should be stable under variations of these 
criteria. 

In the Z-|-jets-|--^T channel, lepton \pt\ spectra are 
used to model neutrino \pt\ spectra. There are sev- 
eral sources of the systematic uncertainty associated with 
this modeling as MC is used to obtain corrections to the 
charged lepton \pt\ spectra. Because these corrections 
are small, uncertainties due to MC used to extract them 
enter only at second order. They can be estimated by 
varying the composition of the MC samples used to mea- 
sure them and the reconstruction efficiencies of leptons 
and jets within their uncertainties. The uncertainties in 
the composition of the MC samples should come from an 
independent measurement of the relative VF-|-jets and 
tt-fjets cross sections for different Nj. Note, in sec- 
tion IIX A[ it is demonstrated that these corrections may 
be negligible for Nj > 4 to obtain a prediction at high 
l^T to 20% or better in early data. 

The QCD background to signal events with one or 
more fake leptons and cross-feeds among F-|-jets pro- 
cesses are not considered in this paper. These back- 
grounds need to be measured and accounted for in the 
I^T distributions and their predictions [2^ . A large new 



physics contamination to QCD at large Jt, in general, 
may bias the prediction at large I^t and hide a new 
physics contribution to V^+jets. I find that under the 
most optimistic scenarios for new physics cross sections 
such a contamination does not lead to a significant bias. 

Even though the reliance on MC is much reduced in 
this method, MC can be used to validate the method and 
constrain its systematic biases as is done in this paper. 
Nevertheless, a study of control data samples is needed 
to develop, optimize and validate the final algorithm and 
to quantify its systematic uncertainties. 

XI. PREDICTIONS WITH SIGNAL 

The algorithm's performance with a new physics con- 
tribution is illustrated in Figure [12] in the .Z-l-jets and 
Z-|-jets+_^T channels for the Nj = 3 and > 4. The inte- 
grated luminosity of the mock data samples in this Fig- 
ure is 200 pb"^ for \/s = 14 TeV. The jet \pt\ thresh- 
olds for Nj is 50 GeV. New physics contributions in the 
Figure are similar to those from mSUGRA benchmarks 
LM4 and LMl ^ for Z-|-jets and ?-|-jets-|--^T, respec- 
tively. The plots show SM backgrounds with new physics 
contributions (dashed) and their I^t predictions (solid) 
from QCD templates. The dot-dashed lines represent SM 
backgrounds only to ease comparisons. 

New physics events tend to have large Jt and I^t- It 
is seen that the addition of a signal contribution with 
large Jt does not bias the prediction significantly at high 
I^T An excess of signal events above the background 
prediction stands out clearly in both channels. Since in 
l+}ets+I^T the neutrino spectrum is modeled based on 
the charged lepton spectrum in each Nj bin, the method 
works best in this signature when the charged lepton 
spectrum in new physics events is soft compared to the 
\pt\ spectrum produced by new weakly interacting par- 
ticles |24| . 

XII. CONCLUSION 

I have presented a new method to predict SM back- 
grounds at high ]^T and a large number of jets, Nj, 
within a context of a search for new phenomena in final 
states consistent with SM F-l-jets: Z+jets, 7-l-jets and 
M^-|-jets. The artificial I^t in each F-l-jets event is 
modeled in-situ using multi-jet QCD events with a con- 
figuration of jets similar to that in the y+jets event. 
The genuine I^t contribution from neutrinos in the 
Z-|-jets-|-^T channel, dominated by T4^-|-jets and tt-|-jets, 
is modeled based on the charged lepton spectrum in each 
Nj bin. 

The method performs reasonably well in robustness 
tests. I have identified mechanisms by which it may 
become biased, discussed systematic uncertainties in its 
background predictions and procedures to estimate them. 
A new physics contamination of the QCD sample does 
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FIG. 12: Observed (dashed) and predicted (solid) SM ^+jets (top) and Hjets+^r (bottom) for Nj = 3 (left) and N,j > 
4 (right) with new physics contributions from mSUGRA benchmarks 17]. The dot-dashed lines highlight the SM contributions. 
The \pt\ threshold for Nj is 50 GeV. The plots correspond to 200 pb~^ at = 14 TeV. 



not lead to a significant bias. The method has discrimi- 
nating po"wer to reveal a ne"w physics contribution at high 
l^T and a large number of jets. It can be applied to data 
■with minimal recourse to MC simulation in early LHC 



running when robust data-driven SM background pre- 
dictions play a key role in searches for ne"w phenomena. 
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